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Charge screening effect in metallic carbon nanotubes is investigated in a model including the 
one-dimensional long-range Coulomb interaction. It is pointed out that an external charge which is 
being fixed spatially is screened by internal electrons so that the resulting object becomes electrically 
neutral. We found that the screening length is given by about the diameter of a nanotube. 



I. INTRODUCTION 

Recently carbon nanotubes(CNTs)S have attracted 
much attention from various points of view. Especially 
their unique mechanical, electrical, chemical properties 
have stimulated many people's interest in the analysis of 
CNTsoH. They have exceptional strength and stability, 
and they can exhibit either metallic or semiconducting 
properties depending on the diameter and helicityoD. Be- 
cause of their small size, properties of CNTs should be 
governed by the law of quantum mechanics. Therefore it 
is quite important to understand the quantum behavior 
of electrons in CNTs. 

However in low energy region ( < 10 4 [K]), we do not 
have to examine all electrons in the system but low en- 
ergy excitations near the Fermi level. Low energy exci- 
tations at half filling move along the tubule axis because 
the circumference degree of freedom (an excitation in the 
compactified direction) is frozen by a wide energy gap( 
~ 10 4 [K]). Hence this system can be described as a 1+1 
dimensional system. Furthermore in the case of metallic 
CNTs, the system describing small fluctuations around 
the Fermi points is equivalent to the two components 
"massless" fermions in 1+1 dimensions. This allows us 
to analyze the low energy excitations with sufficient ac- 
curacy by means of technique in quantum field theory. 

Quantum mechanical behavior of the "massless" 
fermions in metallic CNTs is governed by many kinds 
of interactions. Among others the Coulomb interaction 
is the most important interaction and drives the system 
into a strongly correlated system. The purpose of this pa- 
per is to investigate how the Coulomb interaction plays 
the role when the system is perturbed. Particularly we 
try to answer the following basic question: "What is hap- 
pening when we put an external charge on a metallic car- 
bon nanotube?" The external charge which we consider 
in this paper is externally fixed spinless particle. This 
is thought to be a charged impurity in a metallic nan- 
otubeEI. To answer the above question we investigate the 
charging energy and charge screening effectotj. 

This paper is composed as follows. A model Hamil- 
tonian for low energy excitations in metallic nanotubes 
is constructed in Sec. II. In Sec. Ill, we rewrite the ki- 



netic Hamiltonian in terms of current operators. Then 
we examine the Hamiltonian including the long-range 
Coulomb interaction and analyze some effects of the 
Coulomb interaction on charging energy and charge 
screening effect in Sec. IV. In Sec.|v| we consider the 



gate-electron interaction to show the effect of the long- 
range Coulomb interaction against an external pertur- 
bation. Conclusion and discussion are given in Sec. VI. 
An another type of the Coulomb potential is considered 
m App.[A|. We give an alternative simple derivation to 
obtain main results using the spin-charge separation in 
AppM 



II. THE HAMILTONIAN 

Low energy excitations in metallic carbon nanotubes 
consist of four independent fermion fields. Their quanta 
are spin up(|) and down (J,) electrons in the K and K' 
Fermi points. We denote these fields as, 



1pL,i 



ie (U) x (K,K'), 



(2.1) 



where i is a label for the four fermions and these fermions 
are expressed by \1/_r- t , ^fx i , ^ k'^ , ^k> respectively. For 
example, the field expresses the spin up electron 

field in the K Fermi point. Hereafter we use i as an ele- 
ment of the set: S = {K^,K i ,K' v K[} = tf,l)x(K,K') 
as in eq.(2.1). Each fermion field consist of two com- 
ponents spinor that we name them ipL,i and ipR,i- This 
two components structure is due to the specific lattice 
structure of the graphite sheeto. The subscripts L and 
R denote left-handed component and right-handed com- 
ponent respectively. The left and right are defined by 
the eigenvalue of the matrix 7 5 (= 03) and 03 is the z- 
component of the Pauli spin matrix. 

Time evolution of these fields is governed by the com- 
plicated Hamiltonian in realistic circumstances. For ex- 
ample- tubes are not strict straight line but generally 
bendlilj and the end of a tube(cap) plight mix the wave- 
functions of different Fermi pointal 1 ]. We restrict our 
attention to the most basic and important interactions: 
kinetic interaction, one-dimensional long-range Coulomb 
interaction and gate-electron interaction with a simple 
boundary condition. The kinetic term is 



P x 
—P x 



*i, (2.2) 



where P x = ihd x is the momentum in the tubule axis 
direction(= x G [0 : L] = D) and v F is the Fermi veloc- 
ity. It is remarkable that the dispersion relation is linear 



2 



due to the special band structure of the 7r electrons in the 
graphite sheet. Therefore we call these fermions "mass- 
less" fermions. 

The one-dimensional long-range Coulomb interaction 

is 



J°(x)J°(y) 
8tt J d yj\x-y\ 2 +d 2 



dy, 



(2.3) 



where e is the electron charge and d in the denomina- 
tor denotes the diameter of a nanotube. J (x) stands 
for the sum of each fcrmion charge density Jf(x)(= 



J°(x)=YJ°(x). 



(2.4) 
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A comment is in order about i ntro ducing the cutoff d in 
the Coulomb potential in eq^^). When we write the 
Coulomb interaction without the cutoff, it has an ultra- 
violet divergence in the limit of x — > y. To avoid the 
divergence, we need to introduce a cutoff. It is appro- 
priate to set it the diameter of a tube because when two 
electrons on a nanotube come from opposite sides of the 
tubule axis direction, they repeal and pass each other. At 
the moment they approach most, there is still a distance 
about the diameter of a nanotube to decrease the energy 
of such evcntrr-. This is an origin of the cutoff. Similar 
kind of cutoflO in the Coulomb potential is derived from 
integrating out the circumference degree of the Coulomb 
interaction. Both potentials have the same behavior at 
long distance scale \x— y\ > O(d). However there is a sig- 
nificant difference at short distance scale \x — y\ < 0(d). 
Some detailed discussions about this potential are given 
in App. [A]. 

We also include the effect of gate voltage on a nan- 
otube. The interaction between the gate voltage and 
electrons is 



H G = eV g - J°(x), 



(2.5) 



where V g is a gate voltage which massless fermions feel. 
It should be noted that the gate voltage is not usually 
equivalent to a voltage of the gate itself because of the 
capacitance of a substrate. We may consider the position 
dependent gate voltage V g (x) as a model for an electri- 
cal contact or as an external perturbation to a nanotube. 
Response of the system against a local gate voltage is in- 
vestigated in Sec.|V|. By comparing the system including 
the long-range Coulomb interaction with the one exclud- 
ing the interaction, we show how the long-range Coulomb 
interaction changes the response to the local gate voltage. 
The microscopic Hamiltonian density for this system 

is 



1~t = "Hf + He + Kg- 



(2.6) 



We shall analyze this Hamiltonian without any approx- 
imation. As for the screening of external particles 
which do nO|t-|have a spin, we may neglect "nonlinear 
interactions"ll3 such as backscattering processes because 
total charge sector(B6) decouples from these interactions. 



III. QUANTIZATION OF ELECTRON 

We have specified the Hamiltonian density which de- 
scribes the low energy excitations in metallic nanotubes. 
In this section we quantize the massless fermion fields and 
rewrite the kinetic-.Hamiltonian in terms of the bosonic 
current operatorst3. The Coulomb interaction and gate- 
electron interaction will be investigated in later sections. 

The energy eigenfunctions of the first quantized kinetic 
Hamiltonian are given by 



hFlpv 



h F ipr, 



4> n {x) = —fj-e 1r "-f 3 



(3.1) 
(3.2) 
(3.3) 



where e„ is the energy eigenvalues and L is the length 
of a metallic nanotube. The energy eigenvalues are 
quantized by a boundary condition. Here we impose 
the periodic boundary condition on the wave functions 
ip n (x + L) = ?p n (x). This boundary condition yields the 
following energy spectra: 



2nhvF 



n = An, 



(3.4) 



where A is an energy scale which depends only on the 
length of a nanotube. In order to evaluate the numerical 
value of this energy we use the formula for the Fermi 
velocity 



v F 



3ja 
~2h 



c 

343' 



(3.5) 



where c is the speed of light, 7(= 2.7 [eV])H is the over- 
lap(hopping) integral and a(— 1.42 [A]) is the nearest- 
neighbor distance between two carbon atoms. For a 
single- wall nanotube with L = 3 [pm], A is about 
1.2 [meV]. 

Each massless fermion field can be decomposed into 
the left and right-handed component field as 



$>i(x, t) EE * L:J (X, t) + y R 4x, t), 



(3.6) 



where t is the time. The left and right-handed component 
fields consist of the energy eigenfunctions and annihila- 
tion operators of each quantum, 



^L,i(x,t) 



J2 <^n{x)e 1 1 
nez 



E 

n£Z 



b l Mx)e+ l - 



(3.7) 
(3.8) 



where a l n , W n are independent fermionic annihilation op- 
erator of left-handed component of i-fermion and right- 
handed component of j-fermion satisfying the anti- 
commutators, 



{aloft} = {bib*} =&8 V 



(3.9) 
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All of the other anticommutators vanish. 

In 1+1 dimensions, it is possible to construct the 
fermion Fock space by acting bosonic creation opera- 
tors which are bi-linear of the fermion operators on the 
vacuum. Because the Coulomb interaction consists of a 
product of the charge density, it is very convenient to 
rewrite the kinetic term Hp{= §T~Lf) using the bosonic 
charge density operators. For this purpose, it is useful to 
introduce the left and right currents as follows: 

J°(x) = J L {x) + J R (x). (3.10) 

The left and right currents for each fermion is defined as 

J L ,i(x) = & Lii (x)V Lli (x), 

J R 4x) = tf RA (x)y R ,i(x). (3.11) 

We expand these currents by the Fourier modes. We 
start with the left sector. The left current can be written 
as follows: 



i£S n£Z 



(3.12) 



where each components are given by 

3l = / . Ji,i> 3l, 



ies 



meZ 



1° = 

JLA — 



!L,i- 



The Fourier components j£ i (current operators) satisfies 
the following commutation relations on the fermion Fock 
space(current algebra). 



rn \f| _ 



Tidij Onra i 



(3.13) 



We proceed to consider the right sector in the same way. 
The current operators are 



Mx) = £j R)i (x) = £-?«7 e+i22?£ > 



(3.14) 



3r = £ 3r,i 



3r,. 



nez 
mez 



Jr.,, 
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They satisfy the bosonic commutation relations 



[Jr., (J 



R,3> 



: JR,: 



(3.15) 



The left charge Qla and right charge Qra are con- 
served separately because the system is invariant under 
the following two independent global transformations of 
fermion fields, -> .,• 'I';,' - -> e i9R ^ Rti . We 

must rewrite the fermion Hamiltonian such that all the 
matrix elements are the same as the original fermion 
Hamiltonian by means of current algebra. It is well- 
known that the following Hamiltonian has the sawe-^na- 
trix element as the original fermion HamiltoniantJ'O 



where the f/(l) charge Qi and the chiral charge Q^^ for 
each massless fermion is defined by the summention and 
subtraction between the left and right charges, 



Qi — Ql,i + Qr,i, Qbs — Qla — Qr,', 



(3.17) 



The physical meaning of the chiral charge is spatial inte- 
gration of the electric current density evF(JL(x) — Jr(x)) 
in a tube. This charge measures the left-right asymmetry 
of the vacuum and is defined by the difference of the left 
charge and the right charge on the vacuum. Note that 
the vacuum energy in the last equation is non- vanishing, 
that is due to the finite size effect of a nanotube. Here 
(O) means the vacuum expectation value of an opera- 
tor O. We have defined the second quantized vacuum 
by filling the negative energy modes, leaving the positive 
energy modes empty. 



IV. LONG RANGE COULOMB INTERACTION 

The one-dimensional long-range Coulomb interaction 
in metallic nanotubes is given by the following potential 
form 



H c = / H c dx = — 
Jd 8tt 



J°(x)J°(y) 



y\ 



d 2 



dxdy, (4.1) 



where d is the diameter of a nanotube and this is the 
traces of the pseudo one-dimensional nature of carbon 
nanotubes. Total Hamiltonian of this system consists 
of the kinetic Hamiltonian and the Coulomb interaction. 
The kinetic Hamiltonian is written in terms of the cur- 
rent operators as is shown in the previous section. The 
Coulomb interaction can also be rewritten in terms of 
the current operators because total charge density in the 
Coulomb interaction is a sum of the left and right current 
operators, 



J°(x) = £0^+^)7 
nez 



L 



(4.2) 



Hence, by repl acing J°(x) by eq.(fl.2]) in the Coulomb 
interactional ), we obtain 



+rz £ v[n) ((iE)t + j ^ (j£ + (iS)t) ■ (43) 

n>0 

We have introduced the Fourier components of the 
Coulomb potential, 



H,, 



a£ 

+£ 

n>0 



{{Qi 



i 

12 



3h,i + 3r,i) 



. . „ f n , cos(2na;) 
V(n) = 2 dx 

Jo , , > 



(§)' 



(4.4) 



(3.16) where R is the circumference ( R = %d ) of a nan- 
otube. The Coulomb interaction is given by the sum of 
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zero mode V(0) and nonzero modes V(n). The nonzero 
Fourier mode(n ^ 0) can be approximated to the mod- 
ified Bessel function V(n) ~ 2A"o(2n-^) quite well if 
£ «1. 

We can include external charges in the theory by re- 
placing the total charge density J° with the sum of the 
internal charge J° and external c-number charge density 
J° in the Coulomb interaction, 



H c = 



(J (x) + J° x (x))(J°(y) + J» x (y)) 

V\x-y\ 2 + d 2 



dxdy. 



It is straight forward to derive the corresponding expres- 
sion in terms of current operators 



H n = 



8ttL 



(4.5) 



n>0 



where we introduce f3 n = T-rV(n) and define the Fourier 
components of the external charge as follows: 



n6Z 



(4.6) 



Total charge of the external particle is given by the spa- 
tial integration of the c-number external charge density 
which is referred to Q ex , 



J ex { x )dx — j ex 



(4.7) 



We finally obtain the long-range Coulomb interaction 
in terms of the current operators which modifies the 
energy spectra of the original kinetic Hamiltonian that 
we discussed in the previous section. Before we ana- 
lyze the total Hamiltonian quantum mechanically, it is 
important to know the order of magnitude of the long- 
range Coulomb interaction. Typical energy scale of the 
Coulomb interaction(/3„) of this system is strong as com- 
pared with the energy scale of the kinetic interaction. 
This is shown by the ratio of the energy scale of the 
Coulomb interaction to the energy separation of the ki- 
netic interaction, 



Pn Oi C R 

-r = K {2n-), 

A 7T VF Li 



(4- 



where a(= , e fc ~ -As) is the fine-structure constant. 

V Anne 137 ' 

We used the approximation formula V(n) ~ 2Ko(2n^) 
in eq.(4.8). The value of the ratio is more than one for 
n < 0(js). Therefore this system is thought to be a 
strongly correlated system. So, we should first estimate 
the Coulomb energy scale of the physical quantities in 
this system. For this purpose, we treat the Coulomb in- 
teraction roughly in the next subsection. Detailed anal- 
ysis of the system will be given in latter subsection. 



A. Order estimation 

When an externally fixed particle with total charge 
Q ex is put on a metallic nanotube, the charge density of 
the particle may be modeled by the delta-function, 



J ex( X ) = QexS(x - Xq), 



(4.9) 



where we put the particle at xo(G [0 : L]). It can be 
thought that this external charge distribution is a model 
for a charged spinless impurity. 

It is important to estimate the order of magnitude of 
the energy which we need to put a fixed external charge 
on a metallic nanotube. We simply neglect the interac- 
tion between internal electrons and the external charge, 
so we set (J°(x)} = 0. In this case, we get 



(Hr 



J°,x{x)J°ex{y) 



dxdy 



1 



-,Qi 



8tt JJ d y/\ x - y\ 2 + d 2 " 8tt d 

Estimation of this value for a metallic nanotube with d — 
1.4 [nm] gives (He) ~ 515 -Q\ m [meV]. As is discussed in 
the later section, this energy is strongly modified by the 
Coulomb interaction between internal electrons and the 
external charge. The internal electrons move or rearrange 
in order to decrease the energy of the system. This means 
that the charge screening effect occurs. 

On the other hand, in an experiment concerning trans- 
port on nanotubes, an electron flows from the electric 
contact and after some period it spreads uniformly in a 
nanotube. Then we have 



1 

L' 



(4.11) 



In this case, the Coulomb self energy is 



(H 



e 2 1 



In 



Q 



2 

ex 



mi ^ 

AtiL \ d J 



) 2 . 

* ex 



(4.12) 



For a metallic nanotube with d = 1.4 [nm] and L = 
3 [^ra], we obtain (Hg ) ~ ^-Q\ x [meV]. It should be no- 
ticed that the eq. ( 4.12] ) is the same formula that is usually 
used in the analysis of the Coulomb blockade phenomena 
of a transport experiment of nanotubes. This energy is 
called "charging energy(i? c )" and is not modified by the 
interaction with internal charges due to it's uniform dis- 
tribution, 



4t7Z 



111 



2L 



(4.13) 



An experimental value of the charging energy for a metal- 
lic single-wall CNT with 3.2[/j,m] in length is about E c = 
3.8 [meV]Ej and is very close to the above estimation. 
This is a verification of the cutoff in the potential of the 
Coulomb interaction. 
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B. Screening effect 

We have so far considered the long-range Coulomb in- 
teraction very roughly. In this subsection we solve the 
system quantum mechanically and examine the charge 
screening effect. To investigate this effect, we shall com- 
pute the charging energy and induced charge distribu- 
tion by an external charge. We start with one massless 
fcrmion case. This case corresponds to an approxima- 



tion that we only consider the system with one mass- 
less fermion, for example, \& k 7 ■ Four massless fermions 
case(metallic nanotubes case) can be obtained by straight 
forward extension of the one massless fermion case. 



1. one fermion case 
Total Hamiltonian is given by 



U = A < f (Q)2 + (Q5)2 - i ) f £ ((tfjttf + usfr-g 



n>0 



V 4 12. 

+ E c «Q) + Q ex f + ^ /?„ ((j£)t + j" + (j« )*) (ft + (j«)t + j« ) 



(4.14) 



n>0 



where we omit the subscript i because we consider one 
massless fermion case. The kinetic term and the Coulomb 
term can be combined as follows: 



H = H 



A 
12 



(4.15) 



n>0 



where 
Ho 



~ [{Q? + (Q 5 ) 2 ] + E c ((Q) + Q ex f , (4.16) 

A((j£) t j£ + (jS) t jS) 

& ((jL) f + JS + (&)*) {fl + (JrY + ■ 



(4.17) 

Hq is the zero mode Hamiltonian whose structure is sim- 
ple and is already diagonalized. Hereafter we choose the 
vacuum with vanishing U(l) charge ((Q) = 0) and the 
chiral charge((Q5) = 0) for simplicity. In this case we 
have H = E rj Q\ x . We diagonarize the Hamiltonian 
H n (n 7^ 0) using the Bogoliubov transformation of cur- 
rent operators, 



Jl \ _ ( coshi„ smhtr, 
UrV ) _ \sinh*„ cosht r 



31 , 



(4.18) 



where 



sinh2t„ = cosh2t r 



E n 



(A + (3 n ) , 



E„=\ x 1+^. 



The Bogoliubov transformed currents ( j£, j^) satisfy the 
same bosonic commutation relations as the original cur- 
rent operators eq.( |3.13| ) and eq.(3.15), 



[Jr, (j£) f ] = n6 n 



(4.19) 



After some calculations we diagonarize H n as follows: 



E„ 



((J2) f 
((jS) f + 



+ 7n0e")*) 



(j£ + 7«&) 
-7n(j«)*) + 



A 2 

J ex ) 



- An 
(4.20) 



where 7„ = sinh 2i„ (cosh i„ — sinhi„). It is easy to find 
conditions of the vacuum \vac\;J® x ) in the presence of 
the external charges("i>etci" in the ket denotes vacuum 
state for 'one' massless fermion.), 

Ql+lnSZ) \vac x ;J° ex ) = Q, 

(iS + 7n(i e n J*)Kc 1 ;J e ° x )=0 ) n>0. (4.21) 

We estimate the energy change between the two 
vacua |uaci; 0) (without an external charge) and 
|i>aci; Jg 3 .)(with an external charge) 



SE = (m Cl ;J e °JiJ(j£j|m Cl ;J e 
A 2 



EcQl 



Je 

n>0 n 



E 

n>0 



1 



Ml 
A 



\JexJ Je 



(vaci;0\H(Q)\vaci\0) 



(4.22) 



where H(J® X ) denotes the Hamiltonian with the external 
charge distribution J® x . If we neglect the Coulomb inter- 
action between internal charges and an external charge, 
as in the previous subsection, the charging energy is given 
by 



n>0 



(4.23) 



Comparing this with the exact result(4.22) one can rec- 
ognize the difference between these results. The effect of 
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internal charge redis tribut ion enters in the denominator 
of second term in eq. ( p~22|) . Therefore if the energy sepa- 
ration of the kinetic spectrum is huge compared with the 
Coulomb interaction, 



<1, 



(4.24) 



the effect of rearrangement of the internal electrons is 
very small so that the charging energy is hardly modified 
from the previous order estimation. However for metal- 
lic nanotubes, the situation is completely different as we 
have mentioned at eq.(4 



Next, let us consider the rearrangement of internal 
electrons. This rearrangement can be shown in the ex- 
pectation value of the charge density operator: 



J°(*) = E (0-2)* 



In) Z e 



E 
E 

n>0 




(4.25) 

The expectation value can be calculated using the def- 
inition of t he vacuum in the presence of the external 
charge (fi~2~i"|), 



(J°{x)) l = (vacuJ^l^lvacuJi, 



n>0 



' Jex e 



M \ n J- 



Mil 
A 



n>0 



f 



2|n L 
A 



ex) c 



3ex e 



(4.26) 



For the delta- function charge distribution ( 4. E ), we obtain 

(J°(ar))i =Qex^- _ x \ ycos(^(a;-a;o)). 
„>o 1 + — L L 



(4.27) 



follows: 

H n = A(( J £, 1 ) t Jli + (jS, 1 ) t jl, 1 ) 
+ A ((iE, 2 )ti£, 2 + (j£, 2 )tj£ )2 ) 

+ &n m,x) ] +r R ,i + (fu f +fh + bzr) 

x (Jli + {Jlj + 32,2 + Wb) f + fix) ■ (4-28) 
First we focus on the fermion labeled by 1 . If we regard 
•7*2,2 + UrjV +3 ex as j ex in the previous analysis, we get 
the following expression instead of eq. ( 4.2CQ 



Hn = E n [ (Q2 J + 7n [(j£,3) t + 3R,2 + tiexT] ) 

x(iE,i+7«[tf,2 + CiS l a) t +&]) 

+ ((iS,i) + + 7n[i2,2 + (iS, 2 ) t +ie"x]) 

An 



x (JS,i + 7n [(j£, 2 ) f + iS,2 + OTJ*] ) + n 

+a (0-2,2)^2,2 + 01,2)^,2) 

[(j-2, 2 ) t + iS,2 + (??*)*] |>'2 ia + C?*,^ , 

(4.29) 

The last two lines in the above equation have almost 
same structure of t he one which we have analyzed in the 
previous subsection J4.f 7| ), 

A (0/2,2)^2,2 + ($,2)^,2) 

+Pn% [{31,2^+31,2 + (&,)*] [i2,2 + (Jl, 2 ) f + j£j ■ 

(4.30) 

However in the present case, the Coulomb interaction 
term is modified as follows: 



A 2 

fin —> Pnjjp- 



(4.31) 



We use the same method (Bogoliubov transformation) 
to diagonarize this part of the Hamiltonian with rotation 
angle s n which is different from the previous one(i n ). 



3l,2 
(31,2V 



I coshs r , 
y sinh s„ 



sinh s n 
cosh 8, 



3l,2 
(3I2Y 



,(4.32) 



where 



Pn%2 \f A 2 \ 

smh2s„ = — cosh2i„ = — A + f3 n -=^ I , 

i*n i*n \ ^n J 



2. two fermions case 

Let us proceed to the two fermions case. Here we use 
i = {1, 2} instead of, for example, i = {iff, K{\ for sim- 
plicity. Even in the two fermions case zero mode sec- 
tor of the Hamiltonian is very simple and it gives only 
the charging energy E c Q 2 ex . We should diagonarize the 
nonzero modes of the Hamiltonian. This part is given as 



Thi s tra nsformation gives the result which is similar to 
eg- dHI , 



(02,2^ + ^(^^2,2 + ^ 



{61,2^ + SnJ 

A 2 A 2 



,) (3R,2+5 n {3ex)*)+n 



+^ n £2 p2 \3 ex) 3 



ex ? 



An 
(4.33) 
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where S n = sinh 2s„(cosh s n — sinhs„). We define the 
vacuum for the two fermions case with external charge 
distribution. 

(i2, 2 +<U e ") \vac 2] J° ex }=0, 

■6 n (j? x y)\vac 2 ;J ex )=Q, 



[Jr,2 
Otl + 7r 



[i£ ) 2 + (j«, 2 ) t +ie"])Kc 2 ;Jl) = o, 

[ti n L^+3l2 + ti: x y])\vac 2 ;J° ex )=0. 



From eqs.( Og ),(4.33) and the above definitions of the 
vacuum, the energy change of the two vacuum \vac 2 ]0) 
and \vac 2 ; J®) is evaluated as 




■7 -6 -5 -4 -3-2 -1 1 2 3 4 5 6 7 8 9 10 



x/d 



SE = E c Qi 
= E c Ql 



A 2 A 2 

£ @ n ~P2 ~p2 d ex ) 3e 

n>0 n n 



£ 

n>0 



Z A 



(j e "J*j-. (4.34) 



We see that the effect of 'two' fermions on the charging 
energy is the factor 2 in front of the 2^ in the denomina- 
tor of second term. This means that the charging energy 
decreases more as compared with one fermion case. On 
the other hand, the internal electron charge density dis- 
tribution becomes 



(J°(x)) 



E 

n>0 



9 2/3„ 
Z A 



1 



( (j n y 

2P„ £ yJex) 



+ j'e 



(4.35) 

It can be seen that the effect of the two fermions on the 
induced charge density is also the factor 2 in front of the 



- 



terms. 



3. four fermions case 

Four fermions case can be analyzed by the same pro- 
cedure just like two fermions case. The calculation is 
straight forward but lengthy, so we shall omit them 
here. Instead, we give an alternative simple derivation 
in App.[B|. The formulae for the charging energy and in- 
duced internal charge density are 



SE = E c Q 2 e 



E 

Tt>0 



\ J ex J J e 



' A 



(4.36) 



(A*)>4 = £- 



A 



1 



J<- 



8ft 

n>0 ^ 1 A 

(4.37) 

where the subscript 4 of the vacuum expectation value 
{{J°[x))i) indicates the vacuum of the four fermions case. 
Notice that these equations can be obtained from results 
of one fermion case by replacing the level spacing as fol- 
lows: 



A 



A 



(4.38) 



FIG. 1: Position dependence of the induced charg e de nsity 
in the presence of an external spinless point particle(4.9) that 
has unit charge Q ex — 1. The screening length reaches several 
d(diameter of a nanotube). Th e soli d line is given by the 
formula of the four fe rmion s case(4.40), the dashed line is for 
the one fermion case(4.27). This dashed line corresponds to 
the week coupling case(|-) of the four fermions case. We used 
d = 1.4 [nm] and L = 3 [/jm] in these plots. 



For the del ta-function charge distribution of the external 
charge(4.9), we obtain 



SE 



n>0 A 



I 2 

' ex i 



(4.39) 



(J°(x)) 4 



A 



n>0 



M» L V L 
a v 



(x - x ) 



(4.40) 



The numerical value of the above charging energy is 
about 85 • Q\ x [meV] for a metallic nanotube with d = 
1.4 [nm] and L — 3 [/im]. This value is to be com- 
pared with the previous estimation of the Coulomb self- 
energy of the external charge (He) ~ 515 • Q\ x [meV] 
in eq.( 4.10| ). This indicates that many internal electrons 
are influenced by the external charge and rearranged. To 
confirm this we plot the expectation value of the internal 
charge density in Figjl]. 

Notice that the induced charge density spreads within 
the range of 0(d). This fact c an be recognized by the 
following reason. In eq.(4.40), we can approximately 



regard the summention of n by integral of a = 2n^ 
if the length of a nanotube is huge compared with the 
circumference (5 <C 1). This gives 



Qe 



:K (a) i 



o+ 



1 



■—K (a) R 



cos a 



x 



da. 



Therefore a remaining typical length scale is the diameter 
of a nanotube. Off course, the screening length depends 
on the physical parameter such as the coupling constant 
of the Coulomb interaction^) and the Fermi velocity. 
We also plot eq.(4.27) of one fermion case. The equation 
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0.1 - 

0.2 - 

-0.3 - 

0.4 - 

-0.5 - 

-0.6 - 

-0.7 - 




where 



-7 -6 -5 -4 -3 



io [nm] 



FIG. 2: Diameter dependence of the screening length. The 
scale of the length is given by nm. The solid line is for d = 
1.4[nm] metallic nanotube, the dashed line is for d — 2.1[nm] 
case and the dashed-dotted line is for d = 0.7[nm]. We use 
L = 3 [fira] and Q e x = 1 in these plots. 



corresponds to the week coupling (a — > §■) of the four 
fermions case. The screening length is related to a cutoff 
in the long-range Coulomb interaction, then if we take 
an another type of cutoff, the screening length depends 
on it. 

Let us consider nanotubes which have different diam- 
eter. Fig.|] shows the induced charge distribution( 4.40| ) 
for nanotubes with d = 0.7, 1.4, 2.1[nm], Here we suppose 
that these nanotubes have the same length L = 3[/im]. It 
should be noted that the screening length is proportional 
to the diameter of a nanotube. 



V. GATE-ELECTRON INTERACTION 

In this section we consider the position dependent gate 
voltage and discuss the effect of the Coulomb interaction 
between internal electrons in this system. The gate volt- 
age is expanded in a Fourier series, 



Vg[x) 



nez 



(5.1) 



where g has the dimension of voltage. In this case the 
gate-electron interaction consist of the zero mode and 
nonzero modes, 

E G = egv° g Q + £ £ eg [v n g {j n L ^ + «)*j£,J 

+e 9 [v-jh + (v-nr B J]. (5.2) 

Let us concentrate on the nonzero modes of one 
fermion case without an external charge. The Hamil- 
tonian for the nonzero modes is rewritten as 

'((i2) t + r„«)*)$ + r n t£) 



H n — E n 



gj ) \Jh ' ^n^ g 

+ ((JH) t + r„<) (iS + r n (t£)*)+n 

-2E n T 3 n {v»rv» 



An 
(5.3) 




(5.4) 



The above Hamiltonian gives us definitions of the vac- 
uum, 

(j£ + r„<) \vaa;Vg) = 0, 

(j% + r n (v%)*)\vac 1 ;Vg) = 0, n > 0. (5.5) 

With these definitions we can compute the expectation 
value of the charge density operator. As a simple exam- 
ple, we use a local gate potential which has the Gaussian 
form 



V g (x) ~ g 



-k(x-xo) 2 



(5.6) 



where k decides the size of the local gate potential. Here 
we set k = l/d 2 . This gives the induced charge density 
for four fermions case. 



(J°(x)) 4 = 



A^eg 
* A 



^ 1 + 4^ L 



(2im 
cos I - x Q ) 

(5.7) 

In order to see the effect of the long-range Coulomb inter- 
action on the induced charge density, we plot this func- 
tion with the Coulomb interaction and without it(/3 n = 
0) in Fig.||. The plots in FigJ|| show that the long-range 
Coulomb interaction between internal electrons signifi- 
cantly changes the response of the system to the external 
perturbation. The effect of the Coulomb interaction on 
the induced charge density appears in the denominator 
of the above result. It should be noted that finite width 
of the local gate voltage makes the effective range of n in 
the summention small and high frequency modes of the 
potential ineffective. 



VI. DISCUSSION AND COMMENTS 

In this paper we have analyzed the charge screening 
effect in metallic carbon nanotubes. The significance of 
our present work is put as follows. 



We modeled the Hamiltonian (2. 6) describing the low 
energy excitations in metallic nanotubes and solved the 
system in the presence of the external charge. It was 
found that when we put an external particle on a metallic 
nanotube, the electric charge of the particle is screened by 
internal electrons due to the long-range Coulomb inter- 
action between the particle and internal electrons. The 
Coulomb interaction is strong as comp ared with the en- 
ergy scale of the kinetic Hamiltonian ( |4 . 8| ) . This fact 
makes the quantum mechanical screening complete. 

The screening length is given by about the diameter 
of a na not ube in regard to the long-range Coulomb po- 
tential 4.1). However the length depends on the cutoff 
in the potential of the Coulomb interaction(see App.[A|), 



FIG. 3: Position dependence of the induced charge density 
by the l ocal gate voltage(5.6). The solid line is given by the 
formula(5.7) and the dashed line is given by the same formula 
however without the Coulomb interaction(/3 n = 0). We set 
2|2 = 1 for a metallic nanotube with d = 1.4[nm] and L = 
3[pm]. 



therefore it is important to understand the short distance 
behavior of the potential to find more accurate value of 
the screening length. This is a nontrivial problem be- 
cause the short distance corresponds to the high energy 
region, hence we necessarily examine other bands in ad- 
dition to the linear bands. Anyway the effective range of 
screening(about the diameter of a nanotube or less than 
that) is very small compared with the length of a nan- 
otube, the end of a nanotube(cap) can be thought to be 
ineffective to the screening phenomena. 

The formula for the induced charge d ensity in the pres- 
ence of a point particle is given by eq.(4.40). The sum- 
mention of n in this equation converges due to eq.( ^.8[ ). 
Because of that an extra cutoff of n is not necessaryfel. fn 
substance the summention converges up to n ~ 0{h). 
However, the high frequency modes(n > 0{h)) of the 
long-range Coulomb interaction might be influenced by 
the other bands which do not belong to the 'massless' 
dispersion bands. We would like to make a quantitative 
analysis of these massive bands in a future report. 
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APPENDIX A 

We may use an another kind of the long-range Coulomb 
potential which is derived from integrating out the cir- 
cumference degree of the Coulomb interaction with a cut- 



FIG. 4: Position dependence of the Coulomb potentials. The 
solid line is given by the formula (A2) and the dashed line 
is the potential that we used in the text. Here we take d — 
1.4[nm] and a z = 1.4 [A]. 



off a 2 @. The interaction is given by, 

D 2 



H c = 



8- 



J°(x)V(x-y)J°(y)dxdy, (Al) 



D 



with the potential, 



V(x) 



2 -K 



7. 2 7T 



d 2 



.(A2) 



K (z) is the complete elliptic integral of the first kind and 
the cutoff a z (~ a) denotes the average distance between 
a 2p z electron and the nucleus. It should be noted that 
a new length scale a z in addition to £ (length of a nan- 
otube) and d(diameter) is coming out. This scale modi- 
fies the behavior of the potential at short distance(x — > 0) 
as is shown in Fig.^. 

Hence Fourier components of the potential, 



V{n) 



2 

7T .lu 



2K Q 2n 



R 



sm x 



dx (A3) 



are different from the previous one V(n) ~ 2Kp(2n 



We plot the induced charge density using eq.(A3) in 
Fig.|5|. It can be seen that the screening length is about 
one half the length of the diameter. 

It is important to note that the effective range of n 
summention in the formula of the induced charge density 
becomes large compared with the Coulomb interaction 
in the text. Within this plot, the summention converges 
up to n ~ 0(-^-). So, it is not clear if the screening 
effect can be recognized in the framework of the "low 
energy" excitations. In order to answer such question, 
we introduce a regulator to examine if the high frequency 
modes of the potential give a significant contribution to 
the final resultB. Here we take a simple regulator, 



reg(n) 



1 



exp[n — n*] + 1 



(A4) 
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5 [nm] 



FIG. 5: Position dependence of the in duc ed charge density 
by an external spinless point particle(4.9) with Q ex — 1. 
The solid line is given by the formula for the four fermions 
case(4.40), the dashed line is for the one fermion case(4.20). 
The dotted-dashed line shows the function in eq.(A5) that is 
regularized by the function (Kj) with n* — 3000 ~ C(-§). We 
use the Fourier components( |^43[ ) and take a z — 1.4(~ a) [A], 
d = 1.4 [nm] and L = 3 [/im] in these plots. 



where n* ~ O(^). 
charge density, 



We define the following induced 



reg ( n )8& 2 /27m. 

n>0 ' A ^ 



(A5) 

This function is also shown in Fig.0. We see from this 
figure that the induced charge density oscillates. This 
is due to the regulator and does not have any physical 
meaning. What needs to be emphasized at this point is 
that the screening effect arises from the contribution of 
the low energy region n < n*. Thus it is concluded that 
the charge screening can be analyzed in the low energy 
physics. 

Let us consider nanotubes which have different diam- 
eter. Fig.^ shows the induced charge distribution for 
nanotubes with d = 0.7, 1.4, 2.1 [nm] by means of the 
Coulomb potential (A2). Here we suppose that these 
nanotubes have the same length L = l[/im]. It should 
be noted that the screening length is proportional to the 
diameter of a nanotube. However as compared with the 
previous potential case (Fig.^) the diameter dependence 
is rather week. 



APPENDIX B 



Here we derive the final results eq.(4.36) and eq.(4.37), 
using the operator decomposition into spin and charge. 
First we consider the density operator for K Fermi point. 
We define the following operators for the K point 



Jl,c k 
Jl,S k 



Jl,k l , 
Jl,k 1 , 



FIG. 6: Diameter dependence of the screening length. The 
scale of the length is given by nm. The solid line is for d — 
1.4[nm] metallic nanotube, the dashed line is for d — 2.1[nm] 
case and the dashed-dot ted line is for d = 0.7[nm]. We use 
the Fourier components (A3) and take a z = 1.4(~ a) [A] and 
L = 1 [/jm] in these plots. 



where we omit spatial dependence of these operator for 
simplicity and Jl,c k > Jl,s k express the charge and spin 
operators of left handed fermions in the K Fermi point. 
Similarly for the right-handed sector at K Fermi point, 
we define 



Jr,c k — Jr,Ki + Jr,k l , 
Jr,s k = Jr,k^ — Jr,k l - 



(B2) 



The commutation relations between these operators can 
be calculated by using the commutation relation of orig- 
inal operators( Ji,j, Jr.*)- For example, we obtain the 
following commutation relation 



[ Jl,c k (x), Jl,c k (y)} = 2—d x 6(x - y). (B3) 
For K' Fermi point, similarly we define 



Jl,c k , — Jl.k'^ + Jl,k'^ Jl,s k , — Jl,k\ ~ Jh,K'i 



J, 



R,C„ 



J 



R.Kl 



Je 



J 



R,K' 



J, 



R,K\ 



Using above definitions of the density operators, we de- 
fine the symmetric and antisymmetric combination con- 
cerning two Fermi points using the above spin and charge 
operators for left-handed sector as follows: 

Jl,c+ = Jl,c k + Jl,c kI1 Jl,C- — Jl,c k — Jl,c kI1 
Jl,s + = Jl.s k + Jl,s k ,, Jl,S- = Jl,s k — Jl,s k ,- 

Similarly for right-handed sector, 



Jr.c + = Jr.Ck + Jr,c k ,, Jr.c- = Jr.Ck — Jr.c k ,, 
(Bl) Jr,s+ = Jr,s k + Jr,s k ,, Jr,S- = Jr,s k — Jr,s k ,- 



11 



We get the following kinetic Hamiltonian in term of these 
new density operators 



H F = A 
+A 
+A 
+A 



L f 1 

■g j D ■ Jl,c+(x) 2 + Jr, C+ {x) 2 : dx - — 



L 



1 

12 
1 ' 
12 

Jl,S- {x) 2 + Jb„S- 0) 2 : dx - — 
d 11 . 



Jl,C-(x) 2 + Jr,C-(x) 2 : dx 



Jl,s+ (x) 2 + Jr,s+ (x) 2 : dx - 



Let us define the following operator and this is actually 
equivalent to the charge density operator, 



Jc+ — Jl,c + + Jr,C+, 
J°(x) = J c+ (x). 



(B4) 
(B5) 



Therefore we find that the total charge sector is decouple 
from other operators and this is famous spin and charge 
separation in one dimensional systems. The Coulomb 
interaction is written only by the total charge density. 
Hence we combine the kinetic term and the long-range 
Coulomb term and define C+ sector as follows: 



He, s A 



LP 1 

g J d ■ Jl,c + {x) 2 + Jr, c+ (x) 2 - dx - — 

(J c+ (x) + J° ex (x))(Jc + (y) + J a ex (y)) 



yj\x-y\ 2 +d 2 



dxdy. 
(B6) 



This total Hamiltonian is very similar to the one fermion 
case which we have analyzed. The current operators is 
defined by the same way in the previous section 



Jc + (x) = Y,( i jt,c + ) 



n£Z 



1 



,+»- 



(B7) 



Notice that the commutation relation is modified 

i 

[JL,c+(x),J L<c+ (y)] =4—d x 6{x-y), 
and we obtain the current algebra 

tfl =4n5 nm , 



3l,c+> (j™c+) 



(B8) 

(B9) 
(BIO) 



Analysis of the above Hamiltonian can be done like one 
fermion case, we decompose it into zero mode and 
nonzero modes, 



A 



(Bll) 



where 



n>0 



H = -[(Q) 2 + (Q 5 ) 2 ]+E c ((Q)+Qe X f(B12) 



A 

1 



f(i2,C+)^ 3L,C + + (3R,C + )^jR,C + J 

+ 0n({ji, c y+jR,c + +{ft x y) 

x(i£, c+ + 01,c + ) t +^) ) (B13) 



We can get this Hamiltonian by the following replace- 
me nt in the one fermion Hamiltonian(4.16 4.17), see alos 
eg. dOj , 



A^. 

4 



(B14) 



Therefore formulae of the charging energy and induced 
charge density for four fermions case can be obtained 
by the above replacement in the equations which was 
obtained in the analysis of the one fermion case. 
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